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Abstract 

A low density film on a flat surface is described by an expansion involving the first four virial 
coefficients. The first coefficient (alone) yields the Henry's law regime, while the next three correct 
for the effects of interactions. The results permit exploration of the idea of universal adsorption 
behavior, which is compared with experimental data for a number of systems. 
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I. INTRODUCTION 



Nearly half a century ago the first adsorption experiments were carried out that yielded 
behavior interpreted in terms of two-dimensional (2D) theoretical models Q]. The agreement 
with those models was a consequence of the fact that high surface area forms of graphite 
present extended, uniform and fiat assumed in the models. Subsequently, a wide 

variety of techniques have found similar behavior on other surfaces P|. With an abundance 
of such data comes the opportunity to address questions about universality of the behavior. 
If universal relations are appropriate, it facilitates the planning and interpretation of new 
experiments. 

In this paper, we address this question and several others. To be specific, we examine 
the nature of low coverage adsorption, using the technique of the virial expansion. In this 
case, the expansion includes both gas-surface and gas-gas interactions. Fortunately, virial 
coefficients incorporating the interparticle interactions are available, facilitating this work. 
Hence, the principal effort reported here involves using these coefficients along with a simple 
description of the motion perpendicular to the surface. The results include numerical values 
of the pressure threshold for adsorption and the relationship between the isosteric heat of 
adsorption and the well depth of the gas-surface interaction. In addition, the validity of 
the expansion is tested by comparing the predictions with experimental data for diverse 
adsorption systems. 

II. DESCRIPTION OF MODEL 

The simplest model of a monolayer film in coexistence with vapor considers it as a classical 
noninteracting gas in an external field, the interaction with the substrate 1^. Assuming that 
the adsorbate atoms move freely across the surface, the potential energy of adsorption is a 
function of the normal coordinate, z, alone. An accurate representation of many adsorption 
potentials, sufficient to describe the energetics and dynamics in the vicinity of the minimum, 
is: 




(1) 
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where k is the force constant, D is the well depth and the equilibrium position is z = 0. The 
solution for the number density of a noninteracting gas in the potential V{z) follows as: 

n{z) = n{oo) exp[-pV{z)] (2) 

where f3 = l/lksT), T is the temperature and n(oo) is the vapor density as ^ ^ oo, where 
the true potential energy (but notfl)) vanishes. Then, the film coverage on a surface of 
area A is given by the expression j28|: 

N = aJ n{z)dz = n(oo) J exp[-pV{z)]dz (3) 

With the assumed quadratic dependence ofVonz and the equation for an ideal gas n(oo) = 
f3P, the film's two-dimensional density 6 = N/A is given by: 

e = f3P/{^)expi(3D) (4) 

This Arrhenius form is the generic behavior of adsorption isotherms, with an activation 
energy equal to the well depth. Thus, thermodynamic techniques can be employed for 
finding^he pressure-temperature conditions necessary for adsorption and probe the well 
depth 1^ . However, the noninteracting gas model is applicable only at very low densities. 

An improved model should take into account interactions between the adsorbate atoms, 
the effect of which usually predominates over effects of quantum statistics (except at very 



low temperatures for light and weakly interacting gases) . This improvement is 

provided by the virial expansion, a series expansion whose leading term denotes the ideal- 
system results, while subsequent terms provide corrections arising from the interparticle 
interactions. For an isoenergetic, or structureless, substrate for which the motion is well 
confined to the plane z = 0, the problem becomes essentially 2D. The virial expansion of 
the equation of state of a two-dimensional gas can be written in the form |lol | 

PU = 9 + B2De^ + C2D0^ + D2D9^:. (5) 

where 11 is the 2D (spreading) pressure of the monolayer film. B2d,C2D and D2D are the 2D 
second, third and fourth virial coefficients, respectively. In the 2D case, assumed here, the 
second virial coefficient is related to the two-body interaction through the equation 

/■oo 

B2D = —71" / drr[exp(—/3u{r)) — 1] (6) 
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where r is the intermolecular separation and u{r) is the pair potentiaL The van der Waals 
(VDW) theory is sometimes used to find an approximation of the second virial coefficient 

B^^^ ^ 1/2 / rfV + /3/2 / u{r)d^r (7) 





In the case of a Lennard- Jones type of interaction, 

u{r) = Mi^r - i^)'] (8) 

where e and a are the interparticle well depth and hard-core diameter, respectively. Using 
this interaction, the second virial coefficient in the VDW approximation is 

We computed both the "real" B2D and its VDW approximation, as indicated in Fig.l. Note 
that they differ appreciably at any T. However, one can shift the VDW approximation by 
adjusting the parameters a and b in the equation B2D = b — Pa. The adjustment to a and 
b was done by modifying the interparticle well depth to eynw = l-3e, a correction which 
improves the accuracy of the VDW approximation as seen in the figure. In this way, we got 
a good fit to B2D at intermediate T, corresponding to the range of many experiments. At 
very high T, however, both VDW curves converge to a finite value, while the real B2D goes 
to zero. 

The third virial coefficient is related to the two-body interactions in clusters of three 
particles and is defined by the equation (in the 2D approximation): 

1 r'^ 

C2D = -7T / dr2dr^[eM-l3u{ri2)) - 1] [exp(-/5«(ri3)) - 1] [exp(-/5n(r23)) - 1] (10) 
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This expression omits any explicit three-body interactions j2J]. Details of the computation 
and tables of numerical values of B2D and C2D can be found in the original sources 00. 

The equation of state relates the spreading pressure 11 to ^ and T . However, in an 
experiment one usually manipulates P and T. Therefore, we seek a relationship P{6,T). 
This can be accomplished by employing the equilibrium condition between the gas phase 
and the adsorbed film 

fJ-gas = l^film (H) 

where /i is the chemical potential in the respective phase. The chemical potential of the gas 
(assumed spinless) phase is well known to be 

|3^igas = ln{l3P\^) (12) 



where 

^ " (27rmA;T)i/2 

is the thermal wavelength of the particles. A derivation of fifum appears in the appendix, 
together with a discussion of the limitations of this approach: 

Pfifiim = H^X^) -I3D + 2B2D0 + IC2D0' + ^D^De^ (14) 

Using the equilibrium condition [Eq.ll), the relationship between 6, P and T is found 

e = I3P\ exp(/?D) exp(-2S2Z5^ - ^C^dO^ - ^D^dO^) (15) 

A kind of universal version of this equation appears if the thermodynamic quantities are 
replaced with reduced (dimensionless) quantities. Thus P can be written as a function of 
a dimensionless pressure (-P*): P = P*e/cx^. In the same way, D = D*e, ksT = T*e and 
9 = 9*/a^. The virial coefficients can also be reduced as -82^:) = -820/0"^, C*!/? = C2d/(^^ 



D20 = D2d/<7^- a comparison of Eq.l5 and EqA puts A in correspondence with ^2-11/ Pk, 
which will be used instead of A. The force constant (k) is taken to be D/a^, resulting in a 
universal adsorption equation 



r = py^exp(Z}7T*)exp(-252*z,^* - ^C^*^^^ - Id;^0*') (16) 

An important question (particularly useful for planning experiments) we address is this: 
what value of P* produces "significant" adsorption? By "significant", we mean the opposite 
of non- negligible, arbitrarily defined as reduced density 6* > 0.1. Below this density, the 
virial expansion should suffice over much of the relevant range of T, so we employ it to 
determine this pressure threshold. Fig.2 shows the resulting adsorption threshold curves 
for various numbers of virial terms and for three different values of D* . Note that (for 
a specific well depth) the four curves converge at high T, since then the effect of gas- 
gas interactions disappears. At intermediate T, in contrast, the effect of interactions is to 
enhance the adsorption, which begins at a lower reduced pressure. The limitations of the 
virial expansion approach is seen in the difference between the black and red curves. At 
low T, the curves diverge and the C2D correction is large, as is that of D2D (for T* < 0.45). 
Thus, we conclude that higher order corrections cannot improve the convergence of the virial 
series, since the virial coefficients are themselves divergent for T* < 0.45 

0,0,0- This 

divergence is consistent with the 2D critical temperature of different gases {T* ~ 0.5) 



The curves presented in Fig.2 can be applied to treat specific adsorbates and substrates. 
The standard LJ parameters are well known in the literature jof^l. However, when com- 
paring our results with experiments, we assessed effects of the reduction of LJ well depth 
due to substrate screening 
both sets of parameters. 

Adsorption experiments allow one to find the isosteric heat, which is a rough measure of 
the gas-substrate well depth. The relationship between the isosteric heat and well depth is 
not straightforward. The isosteric heat is 

Qst = -i^)N (17) 

The adsorption equation {Eq.l5) allows us to find the P-T relation, from which one obtains 

In dimensionless form 

^st 2 ^ ^ dT* 2^ ' dT*^ ' ^ ' 

where Q*^ = Qst/^ is the reduced isosteric heat. Fig.?) shows the difference between the 
reduced isosteric heat and the well depth, as a function of T in three cases: noninteracting 
gas, interacting gas including only B2D and interacting gas including both B2D and €20- For 
specificity, we consider the case of 0* = 0.1. Let us focus first at low T: the difference between 
Q*^ and D* is small for the noninteracting model (T*), but for the interacting gas the isosteric 
heat and well depth differ more due to the (overall attractive) gas-gas interactions. (Again, 
caution must be applied with the derivatives of virial coefficients, divergent for T < 0.5.) 
For intermediate values of T*, this difference is of the order of the well depth e, which is 
particularly important in cases when t> D. At T* ^ 2, the interacting and noninteracting 
models yield quite different values of Qst- This can be seen easily for the second coefficient, 
computed in the VDW approximation, 

' dT* ~ 5 ^ 

When T — > 00, however, -By£)VK is still finite while B2D 0. Thus, the noninteracting and 
interacting models converge at high T to the limit Q*^ — D* = T* /2. 
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III. COMPARISON WITH EXPERIMENTS 



A relevant assessment of these expressions is possible for adsorption on graphite and Ag 
(111) substrates, since there is a large body of experimental data for these surfaces. The data 
were obtained from electron diffraction studies of the adsorption of various gases. During 
these experiments, gas was adsorbed under quasi-equilibrium conditions and the intensity of 
the diffraction spots from the substrate was measured. The attenuation of these diffraction 
beams provides a measure of the amount of gas adsorbed; in all of these cases, distinct 
"steps" were observed in the resulting equilibrium isobars or isotherms, corresponding to 
adsorption of the first layer. (Subsequent steps corresponding to subsequent layers were also 
observed in most cases, but those are not relevant to this paper.) By studying this behavios 
over a range of pressure or temperature, it is possible to produce an "isostere" curve in the 
P-T plane for any chosen coverage. The coverages chosen in the experimental studies cited 
here were the half-monolayer and monolayer coverages. Because of the limitation posed by 
the mean-free path of electrons, gas pressures higher than about 10~^ mbar were inaccessible 
in these experiments, and therefore the range of data available for comparison is the pressure 
range between about 10~^^ mbar and 10"'^ mbar. 

As sources for our well depth values, we used Refs. for graphite and [3l for Ag(lll). 
Results of our calculations are shown (using the screened LJ parameters from Table 1) in 



Fig A (for grap hite) and Fiq.p (fo r Agfll 



references lie 



1J,11 




ialialial2i 



)). The experimental data were taken from various 
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[. The calculated curves are in overall good 
agreement with the experimental points, taking into account the approximations we made to 
achieve a simple, universal curve. However, we noticed that the calculations are sensitive to 
the gas interaction parameters; the standard LJ parameters yield less agreement than those 
of the screened interaction. The discrepancies between the calculated curves and some data 
points may be due to the fact that our calculations are performed at coverages much lower 
than a monolayer, while experiments usually tabulate the pressure and temperature values 
related to the completion of the monolayer. However, experimentally, the low coverage 
pressure (at the onset of adsorption) does not usually differ appreciably from the monolayer 
completion pressure. As mentioned above, the virial expansion is not convergent for T < 
0.4e; therefore for gases with large e, some experimental points may be in this range and 
thus other theoretical methods (Monte Carlo simulations) should be applied. 
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IV. CONCLUSIONS 



Virial expansion calculations were performed in order to find the pressure-temperature 
threshold of adsorption on various surfaces assumed to be atomically smooth and perfect so 
they conform to the 2D approximation. Employing the equilibrium condition between the 
film and vapor, an equation relating P to 9 and T was found and reduced to a universal 
equation. A comparison between the nonintcracting and interacting models showed that the 
effect of interactions is to enhance the onset of adsorption, which occurs at lower pressure 
than in the nonintcracting case. At high T, there is little difference between the nonin- 
tcracting and interacting curves, while at intermediate T, the difference between the them 
is larger. At very low T, the virial coefficients are divergent. The computed threshold of 
adsorption was tested for various gases adsorbed on graphite and Ag(lll). The calculated 
curves are in good agreement with the experimental. However, some of the experimental 
points lie in a range where the virial expansion diverges, so the comparison is incomplete. 
Finally, we addressed the relationship between the isosteric heat and the well depth. Due to 
interactions, these quantities differ appreciably at low T, while at high T the nonintcracting 
and interacting cases converge to the same values. 

We thank Mary Jo Bojan and Bill Steele for helpful information and the National Science 
Foundation for support of this research. 

V. APPENDIX 

The film's chemical potential can be found from the thermodynamic equation 



where Qiv is the N-particle partition function of the film. The partition function for a 
nonintcracting 2D gas can be factorized into Q2D, associated with motion parallel to the 
surface, and Qz, associated with motion perpendicular to the surface. Then, the partition 
function can be written as 



where A is the area of the film. The molecules' motion perpendicular to the surface depends 
on the gas-surface interaction potential, here taken as V{z) — —D + kz'^/2. Then, the 




(21) 




(22) 
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quantum partition function for the z motion is 

oo 

Qguan ^ ^ exp[-/3(-D + huj{n + 1/2))] 

n=0 

_ exp(/3(£)-;ia;/2)) 

1 — exp(— /9^a;) 
^ exp(/j^5) 

1 - exp(-/3^c^) ^ ^ 

Here, ^a;/2 is the zero point energy, 1 — exp(— /3^a;) comes from excitations to higher levels 

and Ei, — D — huj/2\s the binding energy. Classically, instead 

dz 



g^' = / ^exp[-/?(-D + ^^V2)] 
exp(/3L') 



Thus, the ratio of Q"^ to 



(24) 



_ 2sinh(/3;ia;/2) 

gr" (/3;k^) ^ ^ 

A high T expansion yields i? f» 1 + (/9/ia;)^/24. This expansion works well up to fihuj = 1, 
at which point R = 1.04. At f3?iuj = 1.5, R = 1.1. Since one takes the logarithm to get 
the chemical potential, a 10% difference in Qz is not significant. For the high T regime, 
Phu < 1.5. Therefore the classical approximation can be used for most purposes, and it 
yields: 

P II film = ln(^A2) -PD + \n{P?iw) (26) 

For (3huj of order 1, the last term in the above equation is negligible, and it makes only loga- 
rithmic corrections elsewhere. Thus, the film chemical potential at low density is essentially 

Plifam^H0oX')-PD (27) 

where = N/A is the film coverage at low densities. In deducing this equation the film was 
supposed to be a 2D noninteracting gas. The calculation of the chemical potential, including 
interactions through the fourth virial coefficient, can be found with the 2D Gibbs-Duhem 

equation 



Integrating 



from a very low density (where the gas can be considered noninteracting) to a higher 
density (9), the film chemical potential becomes 

Pi^fum = In(^A') -(3D + 2B2d9 + ^Csd^' + ^D^dO^ (30) 
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TABLE I: Standard and screened LJ parameters. Notice the well depth (e) reduction due to 
screening when adsorbed on graphite and Ag (111). The change of a due to screening is 

not significant, thus is not shown. 



gas 


^standard \'^) 


Gt (\{'\ 
^screened \ ) 


Ag (XC\ 

^screened \ I 




Xe 


221 


214 


185 


4.10 


Kr 


171 


150 


135 


3.60 


Ar 


120 


110 


96.7 


3.40 


H2 


37.0 


27.3 




3.05 


Ne 


35.6 


33.8 




2.75 


N2 


36.4 






3.32 


CO 


49.6 






3.14 


CH4 


148 






3.45 


He 


10.2 






2.56 
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VI. FIGURE CAPTIONS 



1. The reduced second virial coefficient as a function of reduced temperature. The trian- 
gles correspond to an exact calculation (Eq.6), the circles use the VDW approximation (Eq. 
9) and the square curve is obtained by shifting the VDW curve to the right (corresponding 
to an increase in gas well depth: evDW = l-3e). Inset shows low T behavior. 

2. Reduced pressure at threshold for adsorption as a function of reduced temperature 
for three reduced well-depth values D* . The squares correspond to a non- interacting gas; 
the circles take into account interactions through the triangles include both B2D and 
C2D and the plus curves include the fourth virial term. In each case, little or no adsorption 
occurs below the curve. 

3. Difference between the reduced isosteric heat and well depth as a function of reduced 
temperature in three cases at =0.1: noninteracting gas (triangle), interacting gas through 
B2D (square) and interacting gas with B2D and C2D (circle). 

4. Threshold pressure of various gases on graphite compared with various experiments 
discussed in the text. The dashed curves take into account B2D and the full curves include 
both B2D and C2d- 

5. Threshold pressure of various gases on Ag(lll) compared with various experiments 
discussed in the text. The dashed curves take into account B2D and the full curves include 
both B2D and C2D- 
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